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Abstract. We investigate the laws of coarsening of a two-dimensional system of 
Ising spins evolving under single-spin-flip irreversible dynamics at low temperature 
from a disordered initial condition. The irreversibility of the dynamics comes from 
the directedness, or asymmetry, of the influence of the neighbours on the flipping 
spin. We show that the main characteristics of phase ordering at low temperature, 
such as self-similarity of the patterns formed by the growing domains, and the 
related scaling laws obeyed by the observables of interest, which hold for reversible 
dynamics, are still present when the dynamics is directed and irreversible, but with 
different scaling behaviour. In particular the growth of domains, instead of being 
diffusive as is the case when dynamics is reversible, becomes ballistic. Likewise, 
the autocorrelation function and the persistence probability (the probability that 
a given spin keeps its sign up to time t) have still power-law decays but with 
different exponents. 
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The aim of this paper is to understand the laws of coarsening of a two-dimensional 
system of Ising spins on a square lattice evolving under single-spin-flip irreversible 
dynamics at low temperature from a disordered initial condition. The irreversibility 
of the dynamics comes from the directedness, or asymmetry, of the influence of the 
neighbours on the flipping spin. Yet, at any finite temperature and at long times, 
the system reaches a stationary state whose measure is that of the equilibrium two- 
dimensional Ising model. In other words, the weight of configurations is given by 
the Boltzmann-Gibbs distribution associated to the ferromagnetic Hamiltonian (or 
energy) of the Ising model on the square lattice (see (2.1)). 

We will show that the main characteristics of phase ordering at low temperature, 
such as self-similarity of the patterns formed by the growing domains, and the 
related scaling laws obeyed by the observables of interest, which hold for reversible 
dynamics [1, 2, 3], are still present when the dynamics is directed and irreversible [4, 
5, 6, 7, 8], but with different scaling behaviour. In particular the growth of domains, 
instead of being diffusive as is the case when dynamics is reversible, becomes 
ballistic. Likewise, the autocorrelation function [9, 1, 2, 3] and the persistence 
probability [10, 11, 12] (the probability that a given spin keeps its sign up to time 
t) have still power-law decays but with different exponents. 

A first hint of the occurrence of these phenomena is obtained by measuring the 
speed at which a given minority cluster (a circle, a right triangle, or a square) of (-I-) 
spins in a sea of (—) spins disappears. We thus find that, while the time required 
for the disappearance of this cluster scales as its area for reversible dynamics, this 
time scales as the square root of the area for irreversible dynamics. In other words, 
in this very simple situation, dynamics crosses over from diffusive to ballistic, under 
asymmetric irreversible dynamics. 

We then demonstrate that this also holds for the general situation of a system 
relaxing at zero temperature from an initial disordered configuration. To this end we 
study two complementary facets of the dynamics, namely the relaxation of the energy 
on one hand, and the relaxation of the equal-time correlation function on the other 
hand. 

We first investigate the relaxation of the energy, which yields a first dehnition 
of the typical growing length scale in the system, denoted by ^E{t) (see (4.10)). The 
relaxation of the energy is characterized by two regimes: the scaling regime, where 
coarsening takes place, followed by the late-time regime, where finite-size effects are 
important. In the first regime, ■Cb(I) is observed to grow diffusively (^£;(t) ~ t^/^) for 
reversible dynamics, and ballistically i£,E(t) ~ t) as soon as dynamics is irreversible. In 
the late-time regime the system either reaches the ground state or a blocked metastable 
state. Indeed, though at long times a finite system is expected to reach one of the 
two ground states where spins are either all up, or all down, this is true only with a 
hnite probability. As studied in [13, 14, 15, 16, 17], if one draws a great many Ising 
spin systems on a two-dimensional square lattice, quenched from high temperature 
down to zero temperature and let them evolve under Glauber dynamics [18], then 
a fraction of them gets trapped in blocked configurations, which are stripes running 
across the system. A similar situation holds for our model. A quantitative study of 
this phenomenon is postponed to another work. Here we determine the scale of time 
necessary for the energy to reach its ground state value. We also measure the scale 
of time necessary for the system to reach a blocked state. This allows us to confirm 
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the change of regime in the coarsening behaviour of the system as soon as dynamics 
is irreversible: coarsening becomes ballistic instead of being diffusive as is the case for 
reversible dynamics. 

We then investigate the behaviour of the equal-time correlation function C(r, f). 
This gives access to an alternate definition of a growing length for a system undergoing 
phase ordering, denoted by ^c(f) (see section (5.1)). The ballistic regime is confirmed 
for irreversible dynamics, with ^c(f) ~ b These two facets of the dynamics are 
coherently related since and 5c(f) are essentially proportional. One additional 
finding is the presence of a dynamical anisotropy in the growing length scale ^c{i) for 
any value of the irreversibility parameter V in (2.2). This anisotropy manifests itself 
by a discrepancy between the values of this length according to the direction in which 
it is measured on the two-dimensional lattice. 

We finally proceed to the study of the autocorrelation function and of the 
probability of persistence. Whether power-laws should survive or not in the ballistic 
regime is not obvious to predict a priori. It turns out that both quantities exhibit 
power-law decays, with autocorrelation and persistence exponents larger than their 
counterparts for reversible dynamics. We also study the statistics of the mean temporal 
magnetization which gives another viewpoint on persistence and confirms the values 
of the exponents found by the decay of the persistence probability. The very existence 
of power-law behaviours is a confirmation of the existence of self-similar coarsening, 
even in the presence of ballistic irreversible dynamics. 


2. Definition of the dynamical rules of the directed Ising model 

In this section we first give the general definition of the dynamical rules of our model, 
valid for any temperature, we then specialize this definition to the zero-temperature 
situation under study. 


2.1. Expression of the rate function 


We consider a system of Ising spins on a two-dimensional square lattice of linear size 
L, with periodic boundary conditions. Thus, the number of spins N is equal to L^, 
and the number of bonds is equal to 2L^. The energy (or Hamiltonian), for a given 
configuration of spins, reads 

D = J ^ ^ ^ “f (2-1) 

i—1,L j — l,L 


From now on we will set J = I, = I and denote the reduced coupling constant by 
K = 1/T. 

At each instant of time, a spin, denoted by cr, is chosen at random and flipped 
with a rate denoted for short by w{a). We choose the following form of the rate 
function [4, 6, 7], where {cte, (Tat, ctw, erg} are the neighbours of the central spin cr, 
and E, N, W, S stand for East, North, etc.| 


/ \ IAE, , 1 — V , , 


2 '■ '2 
2l-f E 


<Je<Jn + 7' 


^ 2 — 


( 2 . 2 ) 


X The notation N for North should not be confused with the notation for the number of spins. 
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where a fixes the scale of time, 7 = tanh2_ftr and V, named the velocity, is the 
asymmetry (or irreversibility) parameter, which allows to interpolate between the 
symmetric case {V = 0) and the totally asymmetric ones {V = ±1). This expression 
of the rate function satisfies the condition of global balance, that is to say, leads 
to a Gibbsian stationary measure with respect to the Hamiltonian (2.1) even if 
the condition of detailed balance is not satisfied§. It represents one, among many, 
possible expression of a rate function, for the kinetic Ising model on the square lattice, 
possessing this property. This expression is invariant under up-down spin symmetry. 
In (2.2), for y > 0 (resp. H < 0) the couple (N, E) (resp. (S, W)) is more influential 
on the central spin than the other one. 

The two particular cases of symmetric dynamics (V = 0), and completely 
asymmetric dynamics (V = ±1), deserve special attention. First, if E = 0, then 
the rate function 

2 

w{a) = - [1 - ^Cr(CT_E + <Tn + crw + os) + + O-wOs)] 

(2.3) 


satisfies the condition of detailed balance. This form is different from the usual Glauber 
rate function 

Q/ 

w{a) = - [1 - a tanh[X {as + ctn + ow + cs)]], (2.4) 


which can also be written in terms of spin operators as 
a 


;(cr) = 


2 L 


4(1 -I- 7 ^) 


1 7(2 + 7^) ^ 

^ “ 4(1 + ^ 2 )+ <^N + crw+ os) 

a{aEONOw + onowos + owosoe + osoeon) \ ■ (2.5) 


The Glauber rate function is fully symmetric under a permutation of the neighbouring 
spins, or, equivalently, only depends on the variation of the energy due to a flip, 


AS = 2(t((j_e -I- (Tat -I- ow + crs). (2.6) 

In contrast, the rate function (2.3) with V = Q does not possess this property: it 
is not fully symmetric under a permutation of the neighbouring spins and does not 
depend on AS only. In other words the neighbouring spins of the central spin do 
not all play equivalent roles, i.e., cannot be interchanged. This can be demonstrated 
by looking at the values taken by (2.3) according to the values of the local field 
h = (Je + om + o'w + 0's'- if this sum is equal to ±2 then all rates are the 
same. In contrast, if the sum vanishes then the rate function takes two different 
values, q:(I ± 7 ^), according to the configuration of the neighbours {oe, on, ow,os}' 

a{l 7 ^) corresponds to {-I- H-} or {-h-1-}, while a(I — 7 ^) corresponds to all 

other configurations. This is a manifestation of an anisotropy in the dynamics. By 
comparison, for the Glauber case, if the local field vanishes, the rate function takes 
only one value, a/2. 

Let us add a word of caveat on the terminology used in the present work. As said 
above, the dynamics defined by the rate function (2.3) is anisotropic. Yet, for short, 
when speaking of asymmetric rates we will have in mind the case E 7 ^ 0 . 


§ We refer the reader to [4, 6] for a detailed account of this property. 
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Table 1. Properties of the dynamics investigated in the present work. Glauber 
refers to (2.4), P = 0 to (2.3), and 7 ^ 0 to (2.2). 



Glauber 

V = 0 

V 7 0 

symmetry 

full 

partial 

no 

isotropy 

yes 

no 

no 

directedness 

no 

no 

yes 

reversibility 

yes 

yes 

no 


The case with V = \ corresponds to the totally asymmetric dynamics where the 
central spin a is influenced by its East and North neighbours only. Then 

w{a) = ^ [1 - 7 cr((T_E + an) + 7^ ctectn] ■ (2.7) 

A similar expression, involving aw and as, holds for V = —1: 

“ 70'(o'w + CTs) + 7^ CTivcrs] . (2.8) 

A remarkable fact about these expressions is that they are unique, up to the scale 
of time fixed by the choice of a, in the following sense. For rate functions only 
involving NEC (North, East, Central) spins, or SWC (South, West, Central) spins, 

1. e., when only two neighbours, chosen among the four possible ones, are influential 
upon the central spin, imposing the condition of global balance uniquely determines 
the expressions (2.7) or (2.8) [4, 6]. The rate function (2.7) for the totally asymmetric 
NE dynamics was originally given in [19], under a slightly different form, but without 
considerations on its derivation or its unicity. 

Table 1 summarizes the properties of the dynamics investigated in the present 
work. 

As a final remark, let us point out that the rate function (2.2) is a linear 
combination of (2.7) and of (2.8), which makes its definition rather natural. Of course, 
one could define a more symmetrical rate function by taking a linear superposition 
of the expressions corresponding to the four couples NE (2.7), NW, SW (2.8), and 
SE. This choice would however imply the presence of more than one irreversibility 
parameter in the rate function. 

£.3. Rules of zero-temperature dynamics 

Let us first remind that for zero-temperature Glauber dynamics moves such that 
AE > 0 have zero rates, moves such that AE = 0 have rate 1/2 and those leading to 
AE < 0 have rate 1, in units of a (see (2.4) with 7 = 1). 

In the present case the non-vanishing rates at zero temperature are read off 
from (2.2) and correspond to conhgurations pictured in figure 1. The first line 
correspond to a local field h = 0 {AE = 0), the second line to h = —2 {AE < 0), the 
third one to h = —4 {AE < 0). In units of a, and reading from left to right, the rates 
are equal to 1 — V and 14-V if h = 0, to l — V, l — V, l-\-V, l-\-V, it h = —2, and to 

2, if h = —4. Due to spin-reversal symmetry, the rates are unchanged if one reverts 
the signs of all spins in hgure 1. 

There is a notable difference between Glauber dynamics and the dynamics with 
V = 0, namely that for the former all configurations with h = 0 have rates equal to 
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Figure 1. Local configurations with non-vanishing flipping rates for the central 
spin (in red) in the directed Ising model with rate function (2.2) at T = 0. The 
blue segments represent the domain walls. The first line corresponds to a value 
of the local held /i = 0, the second one to h = —2, the third one to h = —4. 
The associated rates are, from left to right, in units of a, first line: 1 — V^ 1 + V 
(Glauber: 1/2, 1/2); second line: 1 — V^ 1 — V^ IH-V, 1 + V (Glauber: 1, 1, 1, 
1); third line: 2 (Glauber: 1). 



Figure 2. Possible moves for a minority (-I-) square for the directed Ising model 
with rate function (2.2) at T = 0. Only the North-East and South-West corners 
move. The big arrow corresponds to a rate equal to 1 -h y (in units of ck), the 
smaller one to 1 — V (see the first line of figure 1). The same holds for a minority 
(—) square. For Glauber dynamics the flipping rates are the same for the four 
corners. 


1 /2, while, for the latter, configurations with h = 0 where the restricted local field 
(te + O’AT = 0 (and therefore aw + o’g = 0) have zero rates. 

The zero-temperature dynamics with V = 1 (resp. V = — 1) is very simple since 
only one single non-zero rate remains, corresponding to the configurations where the 
North and East spins (resp. South and West spins) are both negative. (For V = —1, 
the same holds for South and West spins.) 

Figure 2 depicts the possible moves of a minority (-I-) square with the rate 
function (2.2). Only the North-East and South-West corners can move, as can be 
seen from the first line of figure 1. We thus infer that a stripe oriented parallel to 
the North-East direction will not move, even when V = 0, whereas a stripe oriented 
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Figure 3. Evolution of a stripe oriented perpendicular to the North-East 
direction for the directed Ising model with rate function (2.2) at T = 0. Stripes 
oriented parallel to the North-East direction are blocked while they evolve for 
Glauber dynamics. 


perpendicular to the North-East direction will, as depicted in figure 3. Examples of 
such configurations will be encountered and commented below (see figures 16 and 17). 

2.3. How to compare two dynamics? 

Since in the course of this work we shall often compare the results obtained with the 
rate function (2.2) to those obtained with Glauber rate function (2.4), we now discuss 
the choice of the time-scale parameter a made in our simulations. 

At infinite temperature (7 = 0) the two expressions (2.2) and (2.4) yield rates 
equal to a/2, hence the two time scales of these dynamics are identical. On the other 
hand, in a simulation at finite (or zero) temperature, if, for practical purposes, one 
wants the rates to be less than 1 , one should fix the value of a according to the largest 
rate, simultaneously for both rates. The largest rates are 

Glauber :f(l + 7 )V(l + 7 ^) x2, 

Directed : (^(1 - 1 - 7 )^ = ^ x 4. (2.9) 

In the simulations presented hereafter the choice a = 1/2 has been done for both 
rates. 

3. Lifetime of a minority cluster 

A first approach to the determination of the scales of time occurring in the coarsening 
process is provided by analyzing the lifetime of a (-I-) cluster in a (—) sea, as depicted 
in figures 4-9. 

Figures 4-6 show the fate of a circle under Glauber dynamics, V = 0 dynamics, 
and V = 1 dynamics, respectively. The two first figures show that the scale of time 
for the disappearance of the circle is the same for Glauber and V = 0 dynamics, but 
that an anisotropy is present in the latter case, while in the former the circle shrinks 
isotropically. Figure 6 shows that, when V = 1, the scale of time for the disappearance 
of the circle is much shorter than for the two reversible dynamics of figures 4 and 5, 
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and that the circle is driven to the shape of a right triangle. Figures 7-9 show that 
the final stages of the evolution of a right triangle are the same as for a circle. 

Figures 10 and 11 give the average time tc needed for a cluster to disappear. 
As demonstrated by these figures, there are two scales of time according to the 
value of V. For reversible dynamics, tc scales as the area A of the cluster, while 
as soon as V is positive this time is proportional to '/A. This gives a first hint that 
asymmetric dynamics is ballistic, a feature which will be confirmed in the sequel by 
several methods. We observe on figure 10 that the asymptotic slope of tc against A 
is the same for all clusters, with common value unity (with the choice made here of 
a = 1/2). 

In contrast, A V = 1, the slope of tc against ^/A depends on the shape of the 
cluster. However a simple rescaling yields figure 11, as follows. Consider a (-I-) cluster 
in a sea of (—) spins. It can be inscribed in a right triangle, as depicted in figure 12. 
The tip of the South-West right angle of the triangle turns out to be the South-West 
corner point appearing in the snapshots. Let us denote by b the distance between this 
point and the hypotenuse of the right triangle. Clusters sharing the same value of b 
should shrink equally rapidly. Hence, in order to compare clusters of different shapes 
we plot tc against b = p^fA, where A is the area, and p the scale factor. For instance 
for a right triangle (lower half of a square with sides of length a parallel to the axes), 
A = a^/2, b = 01^/2, thus p = 1. Likewise for the above mentioned square, A = of, 
b = a\/2, hence p = \/2. For a circle of radius a, p = {1 for the diamond, 

p = 3/2, and for the upper right triangle, p = 2. We thus observe on figure 11 that 
the asymptotic slope of tc against b is the same for all clusters, with common value 
very close to 2y/2. More generally we observe that 


Vtc 


V2b 


(3.1) 


a 

for any H 7 ^ 0 and for a arbitrary. 

The time evolution of a shrinking domain under Glauber dynamics has been 
addressed in [ 20 , 21 ] and we refer the reader to these references for an explanation of 
the slope of G against b observed in figure 10.|| For V = 0 dynamics, since the slope 
is unchanged, it is clear that the dynamics of the cluster is dominated by diffusive 
moves, corresponding to the first line in figure 1 . 

We now present a simple argument in order to explain the above observations on 
figure 11 for the case where V 0. We define a speed v (to be distinguished from 
the so-called velocity V) by b = vtc- Let us consider an infinite interface oriented 
in the North-West to South-East direction made of a regular succession of horizontal 
and vertical domain walls. As is well known, one can describe the evolution of this 
interface by the asymmetric simple exclusion process (ASEP) on a line as follows. An 
horizontal domain wall corresponds to a particle, a vertical domain wall to a hole. 
The current of particles is J = {p — q)p{l — p) where p and q are the rates for a 
jump of the particle, respectively to the right or to the left and p is the density of 
particles. Here p = q;(1 -I- V) and q = q;(1 — V), hence J = aV/2 since p = 1/2 in the 
present case. In a lapse of time r such that each particle has advanced by one step, 
corresponding to a displacement of the interface in the South-West direction equal to 
\/2/2, the total current of particles through a bond in the ASEP is equal to Jr = p. 


More recent developments on related subjects can be found in [22, 23, 24, 25, 26]. 
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Hence v = 'J2I{2 t) and finally 


V = 


aV 

Tf' 


(3.2) 


This simple argument, which accounts for (3.1), demonstrates that the bulk of the 
phenomenon comes from the ballistic motion of the interface oriented in the North- 
West to South-East direction. 





Figure 4. Snapshots showing the disappearance of a circle of minority spins 
when using Glauber dynamics (2.4) at T = 0. The system size is L = 300, and 
the radius of the circle a = L/\. Times are t = 0, f*/3, 2t* jS, 0.95f*, where 
t* = A = 7r(300/4)2 Ri 17671. 






Figure 5. Snapshots showing the disappearance of a circle of minority spins for 
reversible rates (2.3) with y = 0 at T = 0. System size and times are as in 
figure 4. 






L 


Figure 6. Snapshots showing the disappearance of a circle of minority spins 
for totally asymmetric rates (2.7) with V = 1 at T = 0. The circle deforms 
quickly into a right triangle. Here t* = 500 which is the approximate value of tc 
extrapolated from figure 11 for \/A = y^(300/4) r: 133. 

To close, let us mention that rescaling the size of the cluster by the relevant 
time scale {\/i for reversible dynamics, or t for irreversible dynamics) yields limiting 
shapes for the shrinking clusters. This has been thoroughly investigated for Glauber 
dynamics [23, 26]. We report our results for the V = Q and V = 1 dynamics in 
figure 13 together with those for Glauber dynamics. The limiting shape of a cluster 
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Figure 7. Snapshots showing the disappearance of a right triangle of minority 
spins when using Glauber dynamics (2.4) at T = 0. The system size is L = 300, 
and the side of the triangle a = L/2. Here t* = A = L^/8 = 11250. 





Figure 8. Snapshots showing the disappearance of a right triangle of minority 
spins for reversible rates (2.3) with V = 0 at T = 0. System size and times are as 
in figure 7. 





L 


Figure 9. Snapshots showing the disappearance of a right triangle of minority 
spins for totally asymmetric rates (2.7) with F = 1 at T = 0. Here t* = 290 which 
is the approximate value of tc extrapolated from figure 11 for \/A = (300/4\/2) ~ 


106. 


for the V = 1 dynamics is a triangle for obvious reasons. More interestingly, the 
limiting shape for the V = 0 dynamics is the same as that for Glauber dynamics 
(after rescaling the latter by a linear factor equal to 2). In figure 14 the theoretical 
prediction [23, 26] is represented by a dashed curve, indistinguishable from the two 
other curves obtained by simulations. The identity between the limiting curves of 
Glauber and V = 0 dynamics stems from the fact that, in the long-time limit, the 
only relevant moves for the evolution of a cluster are those corresponding to the first 
line of figure 1. 

4. Energy at T = 0 

In this section we start the investigation of the scales of time for a finite system 
undergoing phase ordering after a quench from infinite temperature down to zero 
temperature. We address the following questions: what can be said on the relaxation 
of the system to the ground state, or possibly to other configurations of higher 
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Figure 10. Mean lifetime tc of a cluster of minority spins in a sea of majority 
spins as a function of area for Glauber dynamics as well as for the reversible 
dynamics (2.3) with \^ = 0 at T = 0. The data result from averaging over 8000 
independent histories. 



Figure 11. Mean lifetime of a cluster of minority spins in a sea of majority spins 
as a function of the length b (see text and figure 12) for asymmetric dynamics 
with y = 1 at T = 0. The data result from averaging over 8000 independent 
histories. 


energy? What is the influence of the velocity? This study will confirm the first 
set of observations made in the previous section. 








Dynamics of the two-dimensional directed Ising model 


12 



Figure 12. Cluster inscribed in a right triangle and definition of length b. 



Figure 13. Limiting shapes of a cluster for three dynamics. Glauber; black 
circular-like curve, V = O', red lens, V = I', green triangle. 


We monitor the evolution of the system by measuring (minus) the energy of a 
bond, or density of energy, denoted by e{t). When the system reaches the ground state 
the final value of e(t) is equal to one, therefore its convergence in the late-time regime 
to smaller values is the signature of blocked (or long-lived metastable) configurations. 
We start by considering e(t) for the simple case of a one-dimensional chain of spins, 
as a preparation for the study of the two-dimensional system, that we consider next. 

4.I. Linear chain on a ring 

Let the system size be L, which is also the number of spins N in the chain, as well 
as the number of bonds. The energy of the chain, for a given configuration of spins. 
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Figure 14. Parts of the limiting shapes of figure 13 for Glauber and \/ = 0 
dynamics after rescaling the former by a linear factor equal to 2. The dashed blue 
curve is the theoretical prediction obtained in [23, 26]. 


reads 






(4.1) 


thus (minus) the energy per bond reads 


e = 


E 

L' 


(4.2) 


In the ground state (all spins + or all spins —), we have E = —L hence e = 1. More 

generally we have E = —L + 2n_|_where _is the number of unsatisfied bonds (or 

domain walls), because the cost in energy for the creation of a domain wall in the 
system is equal to 2 units, hence 


e = 1 — 


2n-|_ 

L 


(4.3) 


For instance the excited state of lowest energy for a linear chain on a ring corresponds 
to a pair of domain walls, for which we have E = —L + 4, and e = \ — L. 

Equation (4.3) still holds during the temporal evolution of the system. The inverse 

fraction of unsatisfied bonds _ jL, which is the typical distance between the domain 

walls, gives a natural definition of a growing length ^E{t) as 


Thus finally 




L 

n+_ ■ 


2 

^E{t) 


(4.4) 


e{t) = 1 


(4.5) 
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The expressions above generalize easily to the two-dimensional square lattice. The 
energy is given in (2.1). The energy per bond is now given (in absolute value) by 


e = 


E 


(4.6) 


In the ground state, E = —2L^, hence e = 1. More generally we have E = 
—2L^ -\- 2n^_where n_|_is the number of unsatisfied bonds, hence 


e = 1 — 


n+- 

L2 ■ 


(4.7) 


Consider for example a configuration of the system with one (vertical or horizontal) 
stripe. The number of unsatisfied bonds is equal to n_|_= 2L, hence 


e = 1-(horizontal or vertical stripe). 

L 

More generally, more complex structures will be characterized by 


(4.8) 





(4.9) 


where k = 1,2,.... Such structures, as in figures 16 and 17, will be encountered in 
the sequel. 



t 


1.00 

0.98 

0.96 

0.94 

0.92 



1200 


t 



Figure 15. 300 histories for e(t) with (a): V = 0, (b): V = 0.5, (c): F = 1 for 
a system of linear size L = 100. Figure (d) also depicts the V = 1 case, but on 
another scale of time in order to exhibit trajectories jumping from 62 = 0.96 to 
60 = 1. See figure 16 (right panel) for a snapshot of the system just before the 
jump of one of the two time-trajectories from 0.96 to 1 in figure (d). 
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During the temporal evolution of the system, (4.7) still holds, and, as for the 
chain, we define the growing length fEif) as the inverse of the fraction of unsatisfied 
bonds, n-|_/(2L^). We are thus led to the same formula (4.5) as for the linear chain. 


e{f) = 1 - 


2 

^E{t) ’ 


(4.10) 


where 2/^E{t) is the excess density of energy. 


4-3. Existence of blocked and long-lived metastable configurations 

Figure 15 depicts e{t) for 300 histories, for a system size L = 100, with V = 0, 0.5,1. 
From the observation of these hgures one can draw the following conclusions. Firstly, 
the speed of convergence of the trajectories to their limits depend on V, a feature 
in line with the observations made for the disappearance of a minority cluster in 
section 3. Secondly, since the final value of e{t) should be equal to 1 in the ground 
state, we see that a fraction of the histories do not reach the ground state and therefore 
correspond to configurations of higher energy. We can indeed identify several groups 
of time-trajectories. The first group converges to eo = 1, the ground state, the 
other groups to lower values of e: 0.98, 0.96. These plateau values are in accordance 
with (4.9): Ci = 0.98 corresponds to blocked configurations which are either horizontal 
or vertical stripes; 62 = 0.96 corresponds either to blocked configurations which are 
stripes parallel to the velocity (see left panel of figure 16 for a snapshot of such a 
configuration) or to slowly disappearing configurations (see right panel of figure 16 for 
a snapshot of such a conhguration), which are stripes perpendicular to the velocity. 
Note that smaller values of e can actually be reached if the number of histories is 
larger. For instance, for 4800 randomly selected histories, stripes with e = 63 = 0.94 
are encountered, see figure 17 for an example. They correspond to conhgurations with 
two stripes parallel to the velocity. We also observe on figures 15 that the group of 
histories ending in the ground state (in black) is the largest one, followed by those 
ending on ei (in red), then on 62 (in green). A last observation concerns the longer 
scale of time corresponding to the long-lived metastable states mentioned above. This 
can be seen on figure (d) where after a long time two histories which seemingly had 
ended on 62 = 0.96 finally jump to eg = 1. Note that these ‘jumping’ trajectories do 
exist for any value of V. 

We relegate the quantitative study of the blocked configurations to another work. 

4 . 4 . Scales of time for the coarsening regime defined by the growing length ^E{t) 

We now want to determine the scales of time of the (transient) coarsening regime, 
before e{f) (and therefore ^E{t)) reaches a plateau value. Figure 18 depicts ^E{t) 
(see (4.10)) against t for various values of V. For y = 0 we observe a power-law 
growth of fE{t) as before reaching the plateau value, as for the well-studied 
Glauber case [1, 27]. We then observe, for small values of V, a crossover to a linear 
dependence of ^E{t) in V. 

This crossover from diffusive coarsening to ballistic coarsening is best revealed by 
assuming the scaling form for the ratio 

= (4.11) 

with obvious notations. Figure 19 depicts the scaling function hE{x). For x very 
small hE{x) is constant, corresponding to the diffusive regime, then it crosses over 
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Figure 16. Left panel: snapshot of a blocked configuration with e = 62 = 0.96 
and V = 0. Right panel; long-lived configuration with e = 62 = 0.96 and 
V = 1 which eventually converges to the ground state. This is the configuration 
corresponding, in figure 15 (d), to one of the two rapidly evolving time trajectories 
from 0.96 to 1. The system size is L = 100. 



Figure 17. Snapshot of a (rare) blocked configuration with e = 63 = 0.94 and 
V = 0. The system size is L = 100. 


quadratically to the ballistic regime, when x increases. At larger x the scaling function 
becomes linear in x, i.e., ^E{t) grows as Vt and this linear behaviour extends to even 
larger values of the scaling variable x, as demonstrated in figure 19 right panel. 

4 . 5 . Scales of time for the coarsening regime defined by the time to reach a plateau 

We complement the study done above by measuring the characteristic times for the 
trajectories identihed in hgures 15 to reach a plateau value. We restrict our study to 
trajectories reaching the ground state (e = cq = 1 ) or the blocked state with e = ei. 

We specialize first to trajectories reaching the ground state (cq = 1). As 
demonstrated by figure 15 (d) these trajectories fall into two classes: those reaching 
the ground state directly and those reaching the ground state with a delay, because 
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Figure 18. Typical growing length ^E{t) for V = 0, 0.01,..., 1. System sizes 
range from L = 1400 for F = 0 to L = 2000 for F = 1. 



Figure 19. Scaling function (4.11) showing the crossover from diffusive 
coarsening to ballistic coarsening. Same data as in figure 18 for all range of 
values of X = V\/i. 


they first spend a long time in a long-lived configuration^. The great majority of the 
trajectories belong to the first class. Therefore a simple way of defining a typical time 
avoiding the difficulty of discriminating between the two classes consists in recording 
the distribution of all times tq for reaching the ground state and then measure the 
median of this distribution, denoted by tq, for a given maximal time of observation 
tma.x (and a given number of histories). Indeed the occurrence of trajectories which 
eventually converge to cq = 1 at very late times does not affect the value of tq while 
it would affect the value of the average (tq). 

^ See the right panel of figure 16 which depicts one of the two long-lived configurations of figure 15 (d) 
with 62 = 0.96 (L = 100) eventually reaching the ground state. 
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Figure 20. Median tq of the distribution of times tq to reach the ground state, 
for a given maximal time of observation tmax = 200 000. For every system size 
4000 histories are considered. Plots are against for V = 0 (figure (a)) and 
against L for V = 1 (figure (b)). 


This protocol results in figures 20. For V = Q this time is proportional to as 
it should (diffusive regime). For V = 1 this time is proportional to L: the dynamics of 
coarsening is accelerated, it crosses over from diffusive to ballistic. We observe good 
convergence of tq to a definite value as soon as Wistory > 1000. 



Figure 21. Median tq of the distribution of times tq to reach the ground state 
as a function of system size L, for values of V ranging from 0 to 1. The number 
of histories used for each value of V is equal to 4000 and tmax = 200 000. 

We now measure tq for various values of V and L. We find the curves of figure 21. 
We analyze the data by plotting them with the scaling assumption 

To = L^gAVL). (4.12) 

We predict that for a; = VL small, 5 t( 0 ) is finite, hence tq scales as for x large. 
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Figure 22. Same data as in figure 21, with scaling hypothesis (4.12). 


grix) ^ 1/x hence tq scales as VL. The scaling function grix) describing the crossover 
from diffusive to ballistic is depicted in figure 22. 



Figure 23. Median ^ of the distribution of times t\ to reach the blocked state 
ei, for a given maximal time of observation tmax = 200 000. For every system 
size 4000 histories are considered. Plots are against for V = 0 (figure (a)) and 
against L for V = 1 (figure (b)). 


The same analysis can be performed for the time ri to reach the blocked state ei. 
Figure 23 depicts the median of the distribution of times ri, for a given maximal 
time of observation = 200 000. Plots are against for V — 0 and against L for 
V = 1. The same conclusions hold. 

To summarize, as demonstrated in this section, coarsening is diffusive for 
reversible dynamics with V = 0, as for Glauber dynamics, while it becomes ballistic 
for irreversible dynamics with V 0. 
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The aim of this section is to complement the study made above by the investigation 
of the equal-time correlation function. This allows both to confirm the crossover from 
diffusive to ballistic coarsening as soon as E ^ 0 and to demonstrate the existence of 
an anisotropy in the dynamics of our model. 

5.1. Definition of f,c{t) 

Let us consider the equal-time correlation function 

C{r,t) = {(7ofi{t)(7ij{t)), (5.1) 

where In the following we consider this correlation in different directions. 

The directions parallel to the axes, with either i = 0 or j = 0, will be denoted by (1, 0). 
The North-East and North-West directions, where j = i or j = —i, respectively, will 
be denoted by (1,1) and (1,-1). 



Figure 24. Protocol to extract ^c(i) from C(r,t) (V = 0,L = 1400). 

In the coarsening regime, the correlation function (5.1) allows to define an 
alternate growing length, denoted by Cc(t), as follows. We determine the intersections 
of C{r,t) with a constant line C = Cq. The value used in the following is Cq = 0.5, 
i.e., ^c(t) is obtained as the half-height width of C{r, t) (see figure 24), but we checked 
that other values yield the same results (up to a pre-factor) for the growing length 

We measured the three different lengths f,c{t) corresponding respectively to the 
(1,0), (1,1) and (1,-1) directions, both for Glauber and V = 0 dynamics. For 
both dynamics all three lengths increase as as shown in figure 25. Whereas no 
anisotropy shows up for Glauber dynamics, a direction-dependent pre-factor modifies 
the lengths for V = 0 dynamics. A similar anisotropy is encountered when E ^ 0, as 
seen on figure 26. In the following we restrict to the (1,0) direction for the analysis. 
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Figure 25. Three different lengths extracted from the equal-time correlation 
C{r^t) in three different directions for Glauber dynamics (full lines) and V = 0 
dynamics (dashed lines). The data result from averaging over at least 1000 
independent histories. The system size is L = 1400. 



Figure 26. Three different lengths extracted from the equal-time correlation 
C{r,t) in three different directions for V = 1 dynamics. The data result from 
averaging over at least 1000 independent histories. The system size is L = 2000. 


Figure 27 depicts ^c{t) for various values of V. The crossover to the diffusive 
regime can be seen for small values of V. As for the crossover is best revealed 

by assuming the scaling form for the ratio 


^c{t,V) 
Cc(i, 0) 


hciVVi). 


(5.2) 
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This scaling behaviour is demonstrated in the left panel of figure 28, which is, up to a 
proportionality constant, identical to the left panel of figure 19. Likewise, at large x 
the scaling function becomes linear in x, i.e., ^c{t) grows as Vt. This linear behaviour 
extends to larger values of the scaling variable x, as demonstrated in the right panel 
of figure 28. 



t 


Figure 27. Typical growing length ^c{t) tor L = 0,0.01,...,!. System sizes 
range from L = 1400 for V = 0 to L = 2000 for F = 1. At least 1000 histories 
have been used for the averages. 



Figure 28. Left panel: scaling function (5.2) showing the crossover from diffusive 
coarsening to ballistic coarsening. Right panel: same data as in figure 27 for all 
range of values of x = V\/i. 


We finally investigate the dependence of C{r, t) both in space and time. For 
Glauber dynamics and more generally for reversible non-conserving dynamics, it is 
known that the correlation function (5.1) depends on space and time only through the 
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scaling variable r/^c{t) [1, 27]: 

As shown in figure 29 this also holds for = 0 and V = 1. The three scaling functions 
are slightly different. Note that by the very definition of C{r,t) the three curves are 
constrained to intersect at the point (1,1/2). 



Figure 29. Equal-time correlation function C{r,t) as a function of For 

Glauber and \/ = 0 dynamics the system size is L = 1400, whereas for V = 1 
systems of length L = 2000 have been simulated. The data result from averaging 
over at least 1000 independent histories. 


We finally compared £,E{t) to These two lengths, which are two alternate 

representations of the same reality, appear to be essentially proportional to each 
other, as can be seen on figures 18 and 27, with a constant of proportionality slightly 
dependent on V. 

6. Autocorrelation and persistence 

We finally proceed to the study of the autocorrelation function and of the probability 
of persistence. We also study the statistics of the mean temporal magnetization 
which gives another viewpoint on persistence. The investigation of these quantities 
complement the study of the previous sections, which were only concerned by one¬ 
time observables. The autocorrelation function is a two-time quantity while the 
persistence probability or the mean magnetization are multi-time objects, probing 
the entire history of the system. 

Whether power-laws should survive or not in the ballistic regime is not obvious 
to predict a priori. It turns out that both autocorrelation function and probability 
of persistence exhibit power-law decay in this regime, with autocorrelation and 
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Figure 30. Autocorrelation plotted against 5c(^) Glauber dynamics as well 
as for V = 0 and V = 1 dynamics. The slopes in the long-time regime give the 
autocorrelation exponent A. For Glauber and V = 0 dynamics, 80 000 systems 
with L = 1200 have been simulated. For V = 1 a, total of 100 000 systems with 
L = 2000 were simulated. 


persistence exponents larger than their counterparts for reversible dynamics. The 
very existence of power-law behaviours is a confirmation of the existence of self¬ 
similar coarsening, even in the presence of irreversible dynamics (see also figure 33 
in section 7). 

6 .1. Autocorrelation function 

For the two-dimensional Ising model quenched from infinite temperature down to zero 
temperature the autocorrelation 

C{0,t) = {a{Q)a{t)) (6.1) 

has a power-law decay at large times, 

C{0,t) ^ ( 6 . 2 ) 

with autocorrelation exponent A « 1.26 [9, 28]. Figure 30 shows the results of 
simulations for the autocorrelation plotted against £,c{t)- The results obtained for 
Glauber and 1^ = 0 dynamics are compatible with a common value A « 1.26 for 
the autocorrelation exponent. For the dynamics with V = 1 the value found for the 
exponent, A « 1.85, is markedly larger than for reversible dynamics, showing that 
the system decorrelates more rapidly under ballistic coarsening, as can be intuitively 
understood. 
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6.2. Persistence 

For an Ising spin system evolving after a quench from high temperature down to zero 
temperature under non-conserved dynamics, the probability that a given spin did not 
flip up to time t defines the persistence probability [10, 11]. This probability, denoted 
hy p{t), has a power-law decay at large time, with the persistence exponent 0 [10, 11], 

p{t) ~ r^. (6.3) 

An accurate numerical analysis of this exponent for a two-dimensional system of Ising 
spins has recently been performed in [29], which also reviews the existing literature 
on the subject. 

Here our aim is to investigate, first, whether the anisotropy present in the 
dynamical rules of the H = 0 dynamics changes the exponent 9, and, secondly, whether 
under ballistic coarsening the persistence probability keeps a power-law decay, and if 
so with what value of the exponent. 

Figure 31 depicts the persistence probability p(t) against The behaviour for 

the V = 0 dynamics seems to be the same as for Glauber dynamics, but it is difficult 
to conclude on the sole basis of this figure. The common value of the persistence 
exponent has apparent value 9 « 0.205. For H = 1, a much larger value of the slope 
is obtained, yielding 9 « 0.59. 



Figure 31. Probability of persistence plotted against tor Glauber dynamics 

as well as for V = 0 and V = 1 dynamics. The slopes yield estimates for 6z with 
2 = 2 for Glauber and V = 0, whereas 2 = 1 for V = 1. For Glauber and V = 0 
dynamics the data, which have been obtained for systems of linear size L = 1200, 
result from averaging over 200 independent histories. For V = 1 a total of 320 
systems with L = 2000 were simulated. 
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Related to the concept of persistence is that of mean magnetization. For a system 
of Ising spins the mean magnetization Mf in the time interval (0, t) is defined as [30] 

t 

Mt = ^J (6.4) 

0 

where a{t') is the value of a selected spin at time t'. The integral in the right side of 
this definition is equal to the difference — Tf~, where are the occupation times 
of the selected spin i.e., the times that this spin spends respectively in the (+) or (—) 
states up to time t. These quantities therefore encode the statistics of flips that a 
given spin experiences up to time t. The probability density of reads 

/m(to; t) = -^Prob(Mt > m). (6.5) 

dm 

For a persistent spin, i.e., a spin which never flipped up to time t, the value taken by 
Mt is m = ± 1 , depending on whether the initial value of this spin was +1 or — 1 , and 
Prob(Mt = ±1) is just the persistence probability. Hence the density /m(w;<) has 
discrete components at m = ± 1 , equal to p(t)/ 2 . 

It was shown in [30, 31, 32] that, in the long-time limit, the distribution of M* has 
a limiting form /M(m), which is a U-shaped function with a singularity at m = ±1 
simply related to the persistence exponent 0 as"*" 

/m(to) ~ (1 =F (m^±l). (6.6) 

Thus this limiting process provides a stationary definition of persistence, even at finite 
temperature [31, 32]. 

We performed simulations for the three dynamics under study (Glauber, U = 0 
and U = 1 ) with the aim of complementing the study of the persistence probability. 
Figure 32 depict our results. It yields the predictions 9 « 0.20 for Glauber and U = 0 
dynamics, and 9 « 0.59 for U = 1 dynamics. These values are compatible with the 
estimates obtained from the power-law decay of the persistence probability p{t) in 
figure 31. Note that figure 32 is also in favour of a common value for the persistence 
exponent of Glauber and U = 0 dynamics. 

7. Discussion 

In the present work we have investigated the behaviour of an Ising spin system on 
the two-dimensional square lattice undergoing phase ordering at zero temperature 
after a quench from infinite temperature. This system of spins obeys dynamical rules 
defined by the rate function (2.2), where the parameter V controls the directedness 
or asymmetry of the influence of the neighbours on the flipping spin. For U = 0 the 
dynamics is reversible but not isotropic. As soon as U 7 ^ 0 the dynamics becomes 
irreversible. In the two special cases V = ±1 the only influential spins on the central 
spin are, respectively, the North and East spins (V = 1) or the South and West spins 
(^ = - 1 )- 

The most salient outcome of the present study is the influence of irreversibility 
on phase ordering, which manifests itself by a crossover from diffusive coarsening to 
ballistic coarsening as soon as U 7 ^ 0. Note that this change of regime, from diffusive 

+ The distribution of the mean magnetization can also be studied for the sign of a diffusion field, 
see [30, 33, 34]. For Ising spin systems, extensions can be found in [31, 32, 35]. 
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Figure 32. Probability density fM{m) against (1 — m^). The absolute values of 
the slopes give estimates for 1 — 6 where 6 is the persistence exponent. The system 
sizes are L = 4096, and the data result from at least 500 independent histories. 
The times used are t = 10 000 ~ 48) for Glauber and V = {) dynamics and 

t = 4480 (^cW ~ 400 as obtained by extrapolation of the data of figure 27) for 
V = 1 dynamics. 


to ballistic coarsening, does not exist either for the Ising chain [5] or for the spherical 
model [36], which only experience diffusive coarsening. 

This change of regime was demonstrated by several means, namely 

(i) by the melting of a minority cluster, 

(ii) by the behaviour of the growing length 

(iii) by the measurement of the times to reach the ground state or a blocked 
configuration, 

(iv) by the behaviour of the growing length 

(v) by the measurement of the autocorrelation and persistence exponents, 

(vi) and by the investigation of the statistics of the mean magnetization. 

So doing we discovered that the self-similarity of the coarsening process, which 
holds for the usual reversible single-spin flip dynamics, still holds in the presence of 
irreversibility, as testified by the scaling of the equal-time correlation function and by 
the existence of non-trivial autocorrelation or persistence exponents. The well-known 
patterns of coarsening for the usual Glauber dynamics are changed in the presence 
of an anisotropy in the dynamics. This is illustrated in figure 33, which depicts the 
three situations studied in the present work: isotropic diffusive coarsening (Glauber 
dynamics), anisotropic diffusive coarsening (V = 0), and ballistic coarsening {V = 1). 

The presence of metastable or blocked configurations investigated in [13, 14, 15, 
16, 29] is also observed in the present model. We plan to come back to this issue in 
the future. 
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Figure 33. Snapshots of a system of linear size L = 300 for the three 
dynamics, Glauber (isotropic coarsening), V = 0 (anisotropic coarsening), V = 1 
(ballistic coarsening). Times are chosen such that, from left to right, t;cf) = 
4,8,16,32, corresponding to t = 70,280,1120,4680 for the first two cases and 
t = 39, 78,156, 312 for the last one. 


Let US finally mention that the zero-temperature dynamics of the fully asymmetric 
model with y = 1 is similar to that of the Toom model [37] with zero noise and 
with random sequential dynamics. In this context, the ballistic disappearance of 
a right triangle was noted in [38, 39]. The property of generic nonergodicity (i.e., 
two-phase coexistence over a finite region of the phase diagram) observed in the 
Toom model [37, 38, 39] is thus expected to hold as well for the directed Ising model 
considered in the present work (see also [40, 41]). 
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